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The aim of this study was to examine antibiotic resistance (AR) dissemination in coastal 
water, considering the contribution of different sources of fecal contamination. Samples 
were collected in Berlenga, an uninhabited island classified as Natural Reserve and visited 
by tourists for aquatic recreational activities. To achieve our aim, AR in Escherichia coli 
isolates from coastal water was compared to AR in isolates from two sources of fecal 
contamination: human-derived sewage and seagull feces. Isolation of E. coli was done 
on Chromocult agar. Based on genetic typing 414 strains were established. Distribution 
of £ coli phylogenetic groups was similar among isolates of all sources. Resistances to 
streptomycin, tetracycline, cephalothin, and amoxicillin were the most frequent. Higher 
rates of AR were found among seawater and feces isolates, except for last-line antibiotics 
used in human medicine. Multi-resistance rates in isolates from sewage and seagull 
feces (29 and 32%) were lower than in isolates from seawater (39%). Seawater AR 
profiles were similar to those from seagull feces and differed significantly from sewage 
AR profiles. Nucleotide sequences matching resistance genes blajEM, sail, sul2, tet{A), 
and tet{B), were present in isolates of all sources. Genes conferring resistance to 3rd 
generation cephalosporins were detected in seawater (b/acix-M-i and t'/asHV-i2) and 
seagull feces {blacm-2)- Plasmid-mediated determinants of resistance to quinolones were 
found: qnrSI in all sources and qnrBW in seawater and seagull feces. Our results show 
that seawater is a relevant reservoir of AR and that seagulls are an efficient vehicle to 
spread human-associated bacteria and resistance genes. The E. coli resistome recaptured 
from Berlenga coastal water was mainly modulated by seagulls-derived fecal pollution. 
The repertoire of resistance genes covers antibiotics critically important for humans, a 
potential risk for human health. 
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INTRODUCTION 

Fecal contamination in aquatic environments contributes to the 
spread of human pathogens along with the dissemination of 
antibiotic-resistant bacteria. Overexposure to antibiotics is lead- 
ing to increasing levels of resistance in the human (and other 
animals) commensal microbiota (Austin et al., 1999). The resis- 
tomes of fecal bacteria, once in environmental settings, contribute 
with antibiotic resistance (AR) genes to non-resistant indigenous 
microorganisms (Aminov, 2011; Tacao et al., 2014). In aquatic 
systems the cycle may include subsequent transmission of AR 
to human-associated bacteria (Devirgiliis et al., 2011; Figueira 
et al, 201 1). At the same time, enrichment in AR bacteria is pro- 
moted by the presence of antimicrobials or other contaminants 
in the environment. This facts reinforce the need to identif)' the 
sources of antibiotic-resistant bacteria in aquatic environments 
of human usage (Rosewarne et al., 2010; Gomez-Alvarez et al., 
2012). 



Escherichia coli is a natural member of the intestinal micro- 
biota of human and other homeothermic animals and is a fecal 
indicator bacteria of election (see review by Harwood et al, 2014). 
Multi-resistant strains are found in environments such as coastal 
waters. Although E. coli is considered harmless and commensal, 
several pathotypes have been recognized as significant human 
pathogens, comprising strains implicated in gastro-intestinal, uri- 
nary, or respiratory infections (Kaper et al., 2004). E. coli strains 
can be assigned to four main phylogenetic groups: A, Bl, B2 
and D (Clermont et al, 2000). Extraintestinal pathogenic bacte- 
ria usually belong to groups B2 and D, precisely those that include 
strains with a larger diversity of virulence factors. The commen- 
sal strains belong to groups A and Bl (Picard et al., 1999; Johnson 
etal, 2001). 

In recent years, E. coli has been recognized as a major player 
in the dissemination of AR (Henriques et al, 2006a; Zhao and 
Dang, 2012). Resistance genotypes so far characterized include 
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genes for resistance to last-line antibiotics such as 3rd genera- 
tion cephalosporins and fluoroquinolones (Alouache et al., 2012; 
Wellington et al., 2013; Tacao et al., 2014). To some extent, sim- 
ilar resistance phenotypes and genotypes are shared by strains 
exposed to similar environmental pressures (Wicki et al., 2011). 
Considering this fact, some animal hosts may favor the preferen- 
tial development of certain genotypes. 

Seagulls readily utilize food sources of human origin, especially 
garbage, which may contribute to the high levels of AR associated 
to their microbiota (Bonnedahl et al., 2009; Poirel et al., 2012). 
Multi-resistant phenotypes are found in E. coli strains isolated 
from seagull droppings in different geographic areas (Gionechetti 
et al, 2008; Bonnedahl et al, 2009; Doleiska et al, 2009; Poirel 
et al, 2012). Genes encoding ESBLs (extended spectrum beta- 
lactamases), such as blacix-M-i and blacjx-M-i5 (Bonnedahl 
et al, 2009; Dolejska et al, 2009; Simoes et al, 2010), and resis- 
tance to quinolones (like the plasmid-mediated qnrS) (Literak 
et al., 2010; Vredenburg et al, 2014) integrate the repertoire of 
seagull's resistome. 

Freshwater environments are well-recognized reactors for the 
dissemination and evolution of AR (Baquero et al, 2008; Figueira 
et al., 2011; Tacao et al, 2012) but processes occurring in seawa- 
ter are less understood. The presence of AR strains of E. coli in 
coastal water represents a health issue in areas that are used for 
recreational activities. 

Berlenga is an uninhabited island, recently added to the World 
Network of Biosphere Reserves due to its ecological relevance. 
Episodes of fecal pollution have been recorded in its coastal water 
(Araujo et al., 2014) with seagulls confirmed as the main source 
of fecal contamination. Human-derived sewage (from visiting 
tourists) was also suggested as a secondary source (Araujo et al, 
2014). 

Because it is a well-delimited ecosystem, the island constitutes 
a privileged location to study pathways implicated in dissemina- 
tion of AR in natural environments. Resistance to 3rd generation 
cephalosporins or fluoroquinolones, used as front-line drugs for 
treatment of Gram-negative infections (Canton and Coque, 2006; 
Cattoir et al, 2007), deserved special attention. 

In this study the focus was put in coastal water and in the spe- 
cific contribution of the two previously identified sources of fecal 
pollution. To accomplish our aims, E. coli strains were charac- 
terized at the phenotypic level with respect to their resistances 
profiles; at the genotypic level, specific genes were used as indi- 
cators for the resistance to last-line use antibiotics. Phenotypic 
and genotypic data from coastal water strains was compared with 
data obtained from the strains of each of the two sources of fecal 
contamination. Main contributions to the water resistome are 
hypothesized and discussed. 

MATERIALS AND METHODS 
STUDY SITE AND SAMPLE COLLECTION 

Berlenga Island is the main island of the Berlengas archipelago, 
located in the Portuguese continental shelf, 5.7 miles northwest 
of Cape Carvoeiro. The island has a great diversity of habitats, 
marine species, and avifauna, including a large population of 
the yellow-legged guU {Larus [cachinnans] michahellis). Recently, 
Berlengas was added to the UNESCO'S World Biosphere Reserve. 



Although uninhabited, many tourists visit the island in the sum- 
mer mainly for recreational activities such as swimming and 
diving (Araujo et al., 2014). Human-derived sewage is collected 
on a settlement tank and discharged near the coastline of the 
island. 

Samples of seagull feces, human-derived sewage and sea- 
water from the Berlenga beach were collected as described in 
Araujo et al. (2014). Sampling was performed every 2 weeks 
between May and September 2011. Briefly, samples of seagull 
feces, human-derived sewage and seawater were collected in each 
sampling date and consisted in: (i) 5 composite samples of seagull 
feces scattered on the beach; (ii) 250 mL of raw human-derived 
sewage taken from the effluent of the sanitary infrastructures; 
and (iii) 2 L of seawater collected about 30 cm below the surface. 
Samples were collected in sterile containers and kept on ice until 
processing. 

£. COLI STRAINS 

A total of 939 E. coli isolates were retrieved in Chromocult 
Coliform Agar plates (Merck, Germany) from seagull feces (n = 
427), sewage (« = 170) and seawater (« = 342) as described 
by Araujo et al. (2014). Genetic diversity of the E. coli iso- 
lates was inspected by BOX-PCR with primers and conditions 
previously described (Araujo et al., 2014). Fingerprinting anal- 
ysis allowed the selection of 414 non-clonal isolates, which were 
included in the present study. From those, 179 isolates were from 
seagull feces, 69 isolates were from sewage and 166 were from 
seawater. 

DETERMINATION OF £ COLI PHYLOGENETIC GROUPS 

The triplex PGR developed by Clermont et al. (2000), with 
primers for genes chuA and yjaA and for the DNA fragment 
TSPE4-C2, was used to determine the phylogenetic groups of the 
E. coli strains. Template DNA was obtained by suspending 2 bac- 
terial colonies in 100 |xL of sterile water. For each PGR reaction, 
1 (iL of this suspension was used as template. Appropriate positive 
and negative controls were included in the assay. 

ANTIBIOTIC SUSCEPTIBILITY TESTING 

The E. coli isolates were tested for susceptibility to 16 antibi- 
otics by disk diffusion method on Mueller-Hinton Agar (Oxoid, 
Basingstoke, UK) according to the Clinical and Laboratory 
Standards Institute guidelines (CLSI, 2012). The following disks 
(Oxoid, UK) were used: ampicillin (10 [xg), amoxicillin (10(i,g), 
amoxicillin/clavulanic acid (20/10 |xg), piperacillin (100 (i-g), 
piperacUlin/tazobactam (100/10 |xg), cephalothin (30|xg), cef- 
tazidime (30 |ig), cefotaxime (30 |ig), gentamicin (10 (ig), strep- 
tomycin (10|xg), imipenem (10|xg), nalidixic acid (30|ig), 
ciprofloxacin (5 (ig), tetracycline (30 \Lg), chloramphenicol 
(30 |xg) and trimethoprim/sulfamethoxazole (1.25/23.75 (ig). 
E. coli ATCC 25922 was used as quality control. Isolates were 
classified as sensitive, intermediate or resistant according to 
the CLSI recommendations after 24 h incubation at 37° C. To 
determine the percentage of resistant strains to each antibiotic 
and to define resistance phenotypes, intermediate and resis- 
tant isolates were subsequently grouped in a same resistant 
class. 
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ANTIBIOTIC RESISTANCE GENE DETECTION 

Escherichia coli strains displaying resistance and intermediate 
resistance phenotypes were screened by PGR to detect genes 
conferring resistance to ^-lactams (foZaxEM> blasuY^ bla^jx-u^ 
blaiMP, blaviM, bla^pc, blaoxA-is, blacks, bla AmpCAike), tetracy- 
cline [fef(A), fef(B)], quinolones (qnrA, qnrB, qnrS) and sulfon- 
amides (sull, sul2). The genes tested were chosen among the most 
prevalent in clinical and environmental E. coli isolates. The PGR 
reactions were performed in a MyGycler Thermal cycler (Bio-Rad, 
USA). Primers and PGR conditions are presented in Table 1 and 
Table SI (Supporting Information). The reaction mixtures (25 |iL 
total volume) consisted of 6.25 |xL NZYTaq 2x Green Master Mix 
(2.5 mM MgGl2; 200 |iM dNTPs; 0.2U/|xL DNA polymerase) 
(NZYtech, Portugal), 16.25 |iL of ultrapure water, 0.75 jiL of each 
primer (reverse and forward), and 1 |xL of cell suspension pre- 
pared as described above. The presence of foZflAmpC-like genes was 
inspected using a multiplex PGR as described by Dallenne et al. 
(2010). Negative and positive controls were included in each PGR 
experiment. The negative control differed from other reaction 
mixtures by substituting the cell suspension for the same volume 



of dH20; positive control strains are indicated in Table 1. PGR 
products were analyzed by electrophoresis on a 1.5% agarose 
gel and stained with ethidium bromide. All amplicons obtained 
using primers for foZasHV, ^'octx-m. &^«AmpC-like> qnrA, qnrB, and 
qnrS and PGR products obtained using primers for foZfljEM from 
3rd generation cephalosporins-resistant isolates were sequenced 
(accession numbers KM094197 to KM094211). For this PGR 
products were purified with DNA Glean & Goncentrator (Zymo 
Research, USA) following manufacturer's instructions, and used 
as template in the sequencing reactions. Online similarity searches 
were performed with the BLAST software at the National Genter 
for Biotechnology Information website. 

The genomic contexts of both genes were determined by PGR 
using previously described primers and conditions (Saladin et al, 
2002; Eckert et al., 2006; Tacao et al, 2012). Primers targeted 
either fofacTX-M-i or blacux-i and genes previously found in the 
vicinity of these resistance genes (i.e., IS£CPi for foZacTX-M-i and 
blacuY-2\ ORF477 for fo/flcTX-M-i) (Saladin et al., 2002; Eckert 
et al., 2006; Tacao et al., 2012). PGR fragments were sequenced as 
described above. 



Table 1 | PCR primers and conditions for AR genes amplification. 



Target 


Primer sequence (5'-3') 


Amplicon 
size (bp) 


Program^ 


Annealing 
temperature (°C) 


References 


Control 
strains 


blaoxA-48 


OXA_F:TTGGTGGCATCGATTATCGG 
OXA_R: GAGCACTTCTTTTGTGATGGC 


744 


E 


55 


Poire! et al., 2004 


S. xiamenensis IR34 


blajEM 


TEM_F: AAAGATGCTGAAGATCA 
TEM_R: TTTGGTATGGCTTCATTC 


425 


B 


44 


Speldooren et ai., 
1998 


K. pneumoniae 6T 


blasHV 


SHV_F: GCGAAAGCCAGCTGTCGGGC 
SHV_R: GATTGGCGGCGCTGTTATCGC 


304 


B 


62 


Henriques et al., 
2006a 


K. pneumoniae 2s 


blacjx-M 


CTX_F: SCVATGTGCAGYACCAGTAA 
CTX_R: GCTGCCGGTYTTATCVCC 


652 


A 


55 


Lu etal., 2010 


K. pneumoniae Kp40 




IMP_F: GAATAGAGTGGCTTAATTGTC 
IMP_R: GGTTTAAYAAAACAACCACC 


232 


B 


55 


Henriques et al., 
2006b 


K. pneumoniae KP99cl96 




VIM_F: GATGGTGTTTGGTCGCATATCG 
VIM_R: GCCACGTTCCCCGCAGACG 


475 


B 


58 


Henriques et al., 
2006a 


P aeruginosa NTU-39/00 


b/SKPC 


KPC_F: CATTCAAGGGCTTTCTTGCTGC 
KPC_R: ACGACGGCATAGTCATTT 


538 


B 


55 


Dallenne et al., 2010 


K. oxytoca Ko25 


Wsges 


GES_F: AGTCGGCTAGACCGGAAAG 
GES_R: TTTGTCCGTGCTCAGGAT 


399 


D 


57 


Dallenne et al., 2010 


K. pneumoniae 22 K9 


tet{A) 


tetA_F: GCTACATCCTGCCTTC 
tetA_R: GCATAGATCGGAAGAG 


211 


C 


53 


Nawaz et al., 2006 


£ coli M.I. 10. 2 


tet{B) 


tetB_F: TCATTGCCGACCTCAG 
tetB_R: CCAACCATCACCATCC 


391 


C 


53 


Nawaz et al., 2006 


£ coli M.I. 10.1 


qnrA 


qnrA_F: TTCTCACGCCAGGATTTG 
qnrA_R: CCATCCAGATCGGCAAA 


521 


C 


53 


Guillard et al., 2011 


K. pneumoniae Kp 51 


qnrB 


qnrB_F: GGMATHGAAATTCGCCACTG 
qnrB_R: TTYGCBGYYCGCCAGTCG 


261 


C 


53 


Cattoir et al., 2007; 
Guillard et al., 2011 


K. pneumonia Kpl 


qnrS 


qnrS_F: GCAAGTTCATTGAACAGGGT 
qnrS_R: TCTAAACCGTCGAGTTCGGCG 


428 


C 


54 


Cattoir et al., 2007 


K. oxytoca Ko25 


sull 


sul1_F: CTGAACGATATCCAAGGATTYCC 
sul1_R: AAAAATCCCACGGRTC 


239 


C 


50 


Heuer and Smalla, 
2007 


A. media 


sul2 


sul2_F: GCGCTCAAGGCAGATGGCAT 
sul2_R: GCGTTTGATACCGGCACCCG 


293 


C 


60 


Henriques et al., 
2006a 


E. coli A237 



"A, B, C, D, and E represent different PCR programs (see Table SI for details). 
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STATISTICAL ANALYSIS 

AR phenotypes for each isolate were converted into a numeric 
code for each antibiotic: 0 signifying susceptibility, 1 repre- 
senting intermediate resistance and 2 representing resistance. 
Principal Component Analysis (PCA, using a covariance matrix 
model) was used to explore AR patterns, reducing the mul- 
tidimensional data matrix to a bidimensional biplot fit for 
interpretation, as performed in suchlike studies (Parveen et al, 
1997; Su et al, 2012; Pereira et al, 2013). One-way analy- 
ses of variance (ANOVA) were performed on the PCA sam- 
ple scores (principal components 1 and 2 — PCI and PC2) 
to assess significant differences (p < 0.05) in the AR pro- 
files of E. coli isolates among sampling sources or phyloge- 
netic groups. When applicable, post-hoc Tukey tests followed 
ANOVA tests. ANOVAs and PCA were conducted using Minitab 
and CANOCO for Windows version 4.5 (Scientia Ltd., UK), 
respectively. 

RESULTS 

PHYLOGENETIC DIVERSITY OF £. COLI ISOLATES 

The set of 414 E. coli isolates were assigned to the main phylo- 
genetic groups: A, Bl, B2, and D (Figure 1). Overall, group A 
was the most prevalent (45% of the total number of isolates) fol- 
lowed by group Bl (31%), group D (21%) and group B2 (3%). 
Considering the different sampling sources (e.g., seawater, seagull 
feces and human-derived sewage) distribution among phyloge- 
netic groups was rather similar, except for group B2 isolates that 
were only detected in sewage (10% of the total number of isolates 
retrieved from this source) and in seagull feces (4%). 

ANTIBIOTIC RESISTANCE PROFILES 

The susceptibility patterns of the isolates to 16 antibiotics are 
shown in Figure 2. Resistance (or intermediate resistance) was 
detected to all antibiotics tested and 94% of the isolates were 
resistant to at least one of the antibiotics. The most prevalent 
AR was toward streptomycin (83-100%), followed by tetracy- 
cline and cephalothin in seawater (48 and 42% respectively), and 
tetracycline and amoxicillin in the fecal sources (35 and 34% in 
seagull feces and 23% for both antibiotics in sewage) . Isolates were 




Seawater Seagull faeces Sewage 



FIGURE 1 I Distribution of the E. coli phylogenetic groups (A, B1, B2 
and D) among seawater (n = 166), seagull feces (n = 179) and 
human-derived sewage (n = 69). 



more susceptible to imipenem, 3rd generation cephalosporins, 
gentamicin and the combination piperacillin/tazobactam. 

For the majority of the antibiotics tested (e.g., penicillins, 
cephalothin, tetracycline), resistance rates were higher in iso- 
lates from seawater and seagull feces (Figure 2). Resistance to 
imipenem, 3rd generation cephalosporins and ciprofloxacin was 
most common in sewage. Multi-resistance levels (i.e., resistance 
to antibiotics of at least three different classes) were higher among 
isolates recovered from seawater (39%), followed by isolates from 
seagull feces (32%) and sewage (29%). 

PCA analysis confirmed a clear AR gradient, with less resistant 
isolates (isolates resistant to a lower number of antibiotics) con- 
sistently displaying lower PCI scores (i.e., being located at the left 
of the biplot) and more resistant isolates displaying higher PCI 
scores (i.e., being located at the right of the biplot) (Figure 3). 
Although resistant and multi-resistant isolates were found in all 
sampling sources, significant differences in PCI scores were found 
among sources [ANOVA: _F(2, 411) = 3.7, p = 0.025]. This was 
mostly due to a significantly higher PCI score in seawater isolates 
comparatively to sewage isolates (Tukey test, p < 0.05). No sig- 
nificant differences in terms of the AR profile were found across 
phylogenetic groups, either for PCI [ANOVA: F(2. 397) = 1.59, 
p = 0.206] or PC2 [ANOVA: _F(2, 397) = 0.43, p = 0.649]; group 
B2 was excluded from this analysis due to low representativeness. 

A total of 126 different resistance (or intermediate resistance) 
phenotypes were observed but most (n = 86) were represented by 
only one isolate. Resistance to streptomycin was the most com- 
mon phenotype, identified in 90 isolates, followed by resistance 
to streptomycin and tetracycline (30 isolates) and streptomycin 
and cephalothin (27 isolates). Several multi-resistance pheno- 
types were detected. Among these, 10 phenotypes included 8 
antibiotics and 9 phenotypes included 9 antibiotics. Twenty-four 
phenotypes (Table 2) were common to seawater and to at least 
one of the putative pollution sources (i.e., sewage and seagull 
feces). From these, 45% (n = 11) were common to seawater and 



■ Seawater ■ Seagull faeces Sewage 




AMP AMI AMC PRL TZP KF CAZ CTX [MP CN S NA CIP TE C SXT 
Antibiotics 



FIGURE 2 I Total resistance frequencies to 16 antibiotics of £. coli 
strains isolated from seawater (n = 166), seagull feces (n = 179) and 
human-derived sewage (n = 69). AMR ampicillin; AML, amoxicillin; AMC, 
amoxicillin/clavulanic acid; PRL, piperacillin; TZR piperacillin/tazobactam; KF; 
cephalothin; CAZ, ceftazidime; CTX, cefotaxime; IMR imipenem; CN, 
gentamicin; S, streptomycin; NA, nalidixic acid; CIP, ciprofloxacin; TE, 
tetracycline; C, chloramphenicol; SXT, sulfamethoxazole/trimethoprim. 
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FIGURE 3 I PCA biplot of £. coli isolate scores in terms of their AR 
profile. Arrows represent antibiotics: AMR ampicillin; AML, amoxicillin; AMC, 
amoxicillin/clavulanic acid; PRL, piperacillin; TZP piperacillin/tazobactam; KF, 
cephalothin; CAZ, ceftazidime; CTX, cefotaxime; IMP imipenem; CN, 
gentamicin; S, streptomycin; NA, nalidixic acid; CIP ciprofloxacin; TE, 



seagull feces but were absent from human-derived sewage, while 
only 8% (n = 2) of the phenotypes were exclusively detected in 
both seawater and sewage. 

DETECTION OF ANTIBIOTIC RESISTANCE GENES 

The resistance genes elected for this study were detected in 157 
isolates out of 390 displaying a resistance phenotype (Table 3). 
The most commonly detected genes were fokxEM) sull, sul2, 
tet{A), and fef(B), which were prevalent in 6 out of 8 sampling 
moments (data not shown). The blajEM gene (bZflxEM-i accord- 
ing to sequence analysis of representative isolates) was more 
prevalent in isolates from seawater (69% of 68 isolates resistant 
to penicillins) than in isolates from seagull feces (38% of) or 
sewage (24%). Both fef(A) and sul2 were commonly detected in 
seawater (69 and 64%, respectively) and seagull feces (70 and 
60%), but were less frequent in sewage (31 and 43%). Eighty 
six isolates displaying resistance to quinolones were tested for 
the presence of qnr determinants. Among these, qnrS (identi- 
fied as qnrSl by sequence analysis) was the most prevalent, being 
more frequent in seagull feces (17% of 30 resistant isolates). 
This gene was detected in 5 out of 8 sampling dates. Gene qnrB 
(qnrSi 9 according to sequence analysis) was detected in two sam- 
pling dates in isolates from seawater and seaguU feces and qnrA 
was not detected. Gene fcZflsHV-u was detected in one isolate 
from seawater. Genes blacrx-M and fo/flAmpC-like were detected 
in one isolate each, from seawater and seagull feces respectively. 
Sequencing analysis revealed 100% identity with fcfacTX-M-i and 
blacMY-2- The genomic contexts of both genes were determined 
and results revealed genomic contexts identical to the ones previ- 
ously described in clinical isolates (Kang et al., 2006; Tacao et al., 
2012): ISECPl was detected upstream of both genes and ORF477 
was identified downstream blacTx-M-i- The genes encoding 




' ©▼ o Seagull feces 

] ° °» o ^ Seawater 

-1.5- I ° • Sewage 

h 1 1 1 1 1 1 r^ 

-1.0 -0.5 0.0 0.5 1.0 1.5 2.0 2.5 



Principal component 1 (44.1 %) 

tetracycline; C, chloramphenicol; SXT, sulfamethoxazole/trimethoprim. TZP 
CAZ, CTX, IMP have little influence in the ordination, thus respective arrows 
are not perceptible. Scores are depicted according to the number of antibiotic 
classes that isolates are resistant to (0-6; left panel) or source of origin 
(right panel). 



resistance to carbapenems (i.e., fcZfloxA-48, blamp, bla\/iu, feZflcES 
and fc/flKPc) were not detected in any of the E. coli isolates 
studied. 

A total of 27 different AR genotypes were observed. Among 
these, 7 included only 1 resistance gene, 7 included 2 genes and 
8 included 3 genes. Four genotypes included 4 genes: faZflxEM- 
tet{A)-sull-sul2, detected in 4 isolates from seawater and 5 isolates 
from seagull feces; foZaxEM-fef(B)-5uZi-SMZ2, detected in 2 isolates 
from seawater and 2 isolates from sewage; blaj^M-blasHY-ii- 
tet{A)-sul2 and blajEM-blacTx-M-i-tet{A)-qnrSl, each detected 
in 1 seawater isolate. The genotype blajEM.-tet(B)-qnrSl-sull-sul2 
was detected in 1 isolate from seagull feces. Fourteen genotypes 
were shared between seawater and at least 1 of the putative pol- 
lution sources (Table 4). From these, 50% (n = 7) were shared 
between seawater and seagull feces, exclusively. Only 1 genotype 
was common to seawater and sewage, and absent from seagull 
feces. 

DISCUSSION 

The threats for human health represented by fecal pollution in 
aquatic environments depend on whether the source is human 
or non-human (Field and Samadpour, 2007). Dissemination of 
AR is certainly a risk factor. The particular characteristics of the 
Berlenga Island allowed evaluation of the contribution given by 
human-derived sewage and seagull feces, to the AR patterns found 
in E. coli strains that can be grown and isolated from coastal water 
samples. 

High prevalence of resistance to streptomycin, tetracycline, 
cephalothin, and amoxicillin was detected in all sources. Results 
are in agreement with previous studies that demonstrated a wide 
dissemination of resistance to these antibiotics among E. coli 
isolates from Portuguese clinical settings (Freitas et al, 2014), 
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Table 2 | AR phenotypes detected in seawater and in at least one of the putative pollution sources (i.e., seagull feces and sewage). 



Phenotypes Number of isolates 



Seawater Seagull feces Sewage Total (% of the total number of isolates) 



s 


31 


39 


20 


90 (21.7) 


S, TE 


14 


13 


3 


30 (7.2) 


Kf, S 


22 


5 


0 


27 (6.5) 


AML, Kf, S 


6 


4 


1 


11 (2.7) 


IPM, S 


4 


4 


0 


8 (1.9) 


S, NA 


4 


1 


1 


6 (1.4) 


AMR AML, PRL, KR S, TE 


2 


3 


1 


6 (1.4) 


AMR AML, PRL, KR S, TE, STX 


2 


4 


0 


6 (1.4) 


AML, S 


1 


3 


1 


5 (1.2) 


AMR S 


1 


0 


3 


4(1.0) 


S, NA, TE 


2 


1 


1 


4(1.0) 


AMR AML, PRL, S, TE 


2 


1 


1 


4(1.0) 


AMC, S 


1 


2 


0 


3 (0.7) 


S, TE, STX 


2 


1 


0 


3 (0.7) 


AML, KR S, TE 


2 


1 


0 


3 (0.7) 


AMR AML, KR S 


1 


2 


0 


3 (0.7) 


AMR AML, PRL, S 


1 


2 


0 


3 (0.7) 


AMR AML, PRL, S, TE, STX 


2 


1 


0 


3 (0.7) 


AMR AML, AMC, PRL, KR S, TE, C, STX 


1 


1 


1 


3 (0.7) 


AMR AML, PRL, KR S, NA, CIR TE, C, STX 


2 


1 


0 


3 (0.7) 


AMR AML, AMC, KR S 


0 


2 


0 


2 (0.5) 


KR S, TE 


1 


1 


0 


2 (0.5) 


AMR AML, PRL, KR S, TE, C, STX 


1 


1 


0 


2 (0.5) 


AMR AML, AMC, PRL, IPM, S, TE, C, STX 


1 


0 


1 


2 (0.5) 



The number of isolates from each sampling source displaying each phenotype is presented. 

AMF! ampiciilin; AlVIL, amoxicillin: AMC, amoxiciliin/ciavulanic acid; PRL, piperaclilin; KF, cephalothin; IPM, imipenem; S, streptomycin; NA, nalidixic acid; CIP 
ciprofloxacin; TE, tetracycline; C, chloramphenicol: STX, thnmethopnm/sulfamethoxazole. 



Table 3 | Prevalence of resistance genes detected in the £. coli strains according to the isolation source. 



Antibiotics 


Genes 


Seawater 


Seagul 


1 feces 


Sewage 






Resistant 


Positive 


Resistant 


Positive 


Resistant 


Positive 






isolates 


PGR results 


isolates 


PGR results 


isolates 


PGR results 


3rd generation cephalosporins and/or penicillins 


blajEM 


68== 


47 (69%) 


82= 


31 (38%) 


33== 


8 (24%) 




i>/asHV-12 




1 (2%) 




0 




0 




i>/3CTX-M-l 


6" 


1 (17%) 


8" 


0 


8'= 


0 




blacm-2 




0 




1 (13%) 




0 


Quinolones 


qnrBW 


39 


1 (3%) 


30 


2 (7%) 


17 


0 




qnrSI 




2 (5%) 




5 (17%) 




1 (6%) 


Tetracycline 


fef(A) 


70 


48 (69%) 


63 


44 (70%) 


16 


5 (31%) 




fef(B) 




22 (31 %) 




15 (24%) 




6 (38%) 


Sulfonamide + trimethoprim 


sull 


28 


8 (29%) 


32 


8 (25%) 


7 


2 (29%) 




sul2 




18 (64%) 




19 (60%) 




3 (43%) 



^blajEM snd blasnv were inspected In all penicillin-resistant Isolates (n = 183), 22 of which were also resistant to 3rd generation cephalospohns. 
^blacTX-M snd blacMY were Inspected in isolates displaying resistance both to penicillins and 3rd generation cephalospohns (n = 22). 
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Table 4 | AR genotypes detected in seawater and in at least one of 
the putative pollution sources (i.e., seagull feces and sewage). 



Genotype 


Seawater 


Seagull 
feces 


Sewage 


Total 

(%) 


tet{A) 


21 


20 


3 


44 (11.3) 


WajEM. fef(A) 


13 


6 


2 


21 (5.4) 


tet{B) 


7 


8 


2 


17 (4.4) 


b/ajEM. fef(B) 


7 


3 


2 


12 (3.1) 


WajEM. fef(A), sul2 


4 


5 


0 


9 (2.3) 


WajEM. fef(A), sull, sul2 


4 


5 


0 


9 (2.3) 




3 


4 


1 


8 (2.1) 


WajEM, fef(B), sul2 


4 


1 


0 


5(1.3) 


sul2 


1 


3 


0 


4 (1.0) 


b/ajEM, fef(B), sull, sul2 


2 


0 


2 


4 (1.0) 


fef(A), sul2 


1 


2 


0 


3 (0.8) 


fa/aiEM. sul2 


1 


1 


1 


3 (0.8) 


WajEM. tef(A), qnrSI 


1 


2 


0 


3 (0.8) 


b/ajEM. fef(A), fef(B) 


1 


1 


0 


2 (0.5) 



The number of isolates from each sampling source displaying each phenotype is 
presented. 



animals (Radhouani et al., 2009) and aquatic systems (Pereira 
et al., 2013). For antibiotics mainly used in human medicine such 
as 3rd generation cephalosporins, ciprofloxacin and imipenem, 
levels of resistance were low. An uncommon phenotype combin- 
ing resistance to imipenem and streptomycin and susceptibility 
to penicillins was detected in eight strains, sis of which were 
isolated from sewage. Genes encoding carbapenemases were not 
detected in these strains, an indication that non-enzymatic mech- 
anisms, such as the expression of efflux pumps or porin loss 
(Papp-Wallace et al., 2011), may contribute to the phenotype. 
Concerning 3rd generation cephalosporins, the resistance pheno- 
type was explained only in three out of twenty-two resistant iso- 
lates. Novel or uncommon resistance genes may be present in the 
genomes of the remaining strains. On the other hand, mutations 
in the promoter and attenuator sequences of the chromosomal 
ampC gene may also contribute to reduced susceptibility to these 
antibiotics. These mutations have been described in E. coli isolates 
from different sources, including aquatic systems (Mataseje et al, 
2009). 

Levels of AR in coastal water isolates were expected to be low, 
considering the moderate human impact and the limits to sur- 
vival imposed to E. coli by water salinity. However, experimental 
data revealed resistance levels similar to those previously reported 
for polluted freshwater (Ham et al., 2012; Pereira et al., 2013) 
or wastewater (Figueira et al., 2011). Multi- resistance levels were 
higher in coastal water (39%) than in the two sources of fecal 
contamination (32% for seaguU feces and 29% found in human 
sewage). 

Genes frequently hosted by clinical isolates and sharing similar 
genomic contexts (Saladin et al., 2002; Eckert et al, 2006), were 
also hosted by seawater isolates. Examples are ESBL-encoding 
genes such as blacrx-M-i and fo/flsHV-i2 or plasmid-mediated 
quinolone resistance genes (e.g., qnrB19). It is to point that 
dissemination of these genes or related ones is a matter of serious 



concern in a scenario of lack of innovation in what concerns 
antimicrobial therapies. In fact, 3rd generation cephalosporins 
and fluoroquinolones are front-line therapeutics for the treat- 
ment of Gram-negative infections (Canton and Coque, 2006; 
Cattoir et al, 2007). 

AR features of seawater E. coli strains are mainly dependent 
on their animal source (Wicki et al., 2011). At the same time, 
E. coli strains are able to adapt to non-human environments, a 
process that implies changes in their genetic background (van 
Elsas et al, 2011). Adaptation to the more stressful saline con- 
ditions of coastal waters may favor strains with genotypes that are 
more efficient to deal with other stresses. 

The results strongly suggest that the E. coli resistome found in 
coastal water was mainly impacted by fecal pollution from seag- 
ulls. The AR patterns found on seagull feces and seawater shared 
abundant common elements, while differing significantly from 
the resistance patterns found on sewage. Differences were not 
attributable to a different distribution in terms of phylogenetic 
groups. Based on these observations it is possible to hypothesize 
that the enteric environment of seagulls positively selects E. coli 
strains more prone to survive in marine environments. The same 
data allows proposing resistome analysis as a microbial source- 
tracking tool (see Harwood et al., 2014). Such approach has been 
previously tested, using almost exclusively data resulting from 
the analysis of resistance phenotypes (Parveen et al., 1997; Field 
and Samadpour, 2007; Carroll et al, 2009). Data on resistance 
genotypes add valuable information for the identification of the 
contamination source. The resistance gene pool in the analyzed 
sources was quite stable along the period of the study, except 
for the sporadic detection of infrequent genes (data not shown). 
Our recommendation is toward the use of genotypic data only 
in combination with other data. Our data further confirmed the 
conclusion of Araujo et al. (2014) establishing seaguUs as the main 
source of fecal pollution on Berlenga Island. 

This study highlights the risks associated to seagull-derived 
pollution. This is a critical aspect since wild life-derived pollution 
is frequently neglected due to the assumption that host-specificity 
of pathogens (Field and Samadpour, 2007) is a preventive bar- 
rier, that makes difficult if not impossible, the establishment on 
a new host. Indeed, the contribution of seaguU fecal bacteria 
to AR dissemination constitutes an undeniable risk to humans, 
especially because of the close contact between these birds and 
human waste. Previous studies have confirmed seagulls as a reser- 
voir of AR (Gionechetti et al, 2008; Bonnedahl et al, 2009; 
Dolejska et al., 2009). In our study, besides the high level of multi- 
resistance in seagull feces, we retrieved isolates carrying genes 
conferring resistance to antibiotics that are critically important 
for human medicine. For example, the prevalence of isolates car- 
rying the qnrSl gene was higher in seagull feces than in sewage. 
This plasmid-encoded gene was previously identified in isolates 
from seagull feces (Literak et al, 2010) and is frequently found in 
human pathogens (Cattoir et al, 2007). Likewise, we detected the 
presence of qnrB19 in seawater and seagull feces, a gene which was 
previously detected in animal and human isolates (Dolejska et al., 
2013). Moreover, the plasmid-encoded gene blacMY-2> conferring 
resistance to 3rd generation cephalosporins, was detected in seag- 
ull feces and in a genomic context identical to the one found in 
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clinical isolates (Verdet et al., 2009). Poirel et al. (2012) reported a 
high proportion of CMY-2-positive strains among resistant E. coli 
isolated from seagull feces in Miami Beach. So, data obtained 
in different contexts indicate seagulls as important vectors of 
antibiotic resistant bacteria. 

A previous study in the same Island (Radhouani et al., 2009), 
showed that all isolates were susceptible to carbapenems and 3rd 
generation cephalosporins. These antibiotics are mainly used in 
human medicine (Canton and Coque, 2006; Papp-Wallace et al, 
2011) and resistance to them is known to be stUl limited in unpol- 
luted or moderately polluted environments (Henriques et al., 
2012; Tacao et al., 2012). As reported here, resistance to these 
antibiotics was most common in human-derived sewage, but its 
presence on seagull feces isolates suggests an undesirable progress 
in terms of resistance dissemination. Seagulls based on the Island 
fly to urban areas on the coast for feeding. Like suggested for other 
geographic regions (Bonnedahl et al., 2009), human activities 
may considerably influence the AR patterns associated to seagulls 
in remote areas, such as Berlenga. In fact, it is a matter of concern 
that levels of resistance to 3rd generation cephalosporins have 
been consistently increasing in isolates from clinical settings in 
Portugal (ECDC, 2013). For example, gene blacMY-it which was 
here detected in one isolate from seagull feces, has been detected 
in E. coli from clinical samples in Portuguese health care settings 
(Freitas et al, 2014). 

CONCLUSIONS 

This study confirms seawater as a relevant reservoir of E. coli 
strains adapted to marine environments and carrying AR deter- 
minants that are common to clinical strains. It is evidenced that 
the seawater E. coli resistome in Berlenga is mainly modulated 
by seagulls-derived pollution. This fact confirms seagulls as a 
main spread vehicle of AR bacteria to coastal water. Occurrence 
of high levels of multi-resistance and of clinically relevant genes 
underlies the need to carefully consider the human health risks 
subjacent to seagull-derived pollution. The assessment of these 
risks is essential for an effective management of coastal water used 
for recreational activities. Lastly, we believe that this study — along 
with correlated data from previous studies — has succeeded in 
showing that AR could be used as a tool to aid in the identification 
of the fecal sources of contamination. 
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